Correlation dynamics after short-pulse photoassociation 
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Two atoms in an ultracold gas are correlated at short inter-atomic distances due to threshold 
effects where the potential energy of their interaction dominates the kinetic energy. The correla- 
tions manifest themselves in a distinct nodal structure of the density matrix at short inter-atomic 
distances. Pump-probe spectroscopy has recently been suggested [Phys. Rev. Lett. 103, 260401 
(2009)] to probe these pair correlations: A suitably chosen, short photoassociation laser pulse de- 
pletes the ground state pair density within the photoassociation window, creating a non-stationary 
wave packet in the electronic ground state. The dynamics of this non-stationary wave packet is 
monitored by time-delayed probe and ionization pulses. Here, we discuss how the choice of the 
pulse parameters affects experimental feasibility of this pump-probe spectroscopy of two-body cor- 
relations. 

PACS numbers: 32.80.Qk,03.75.Kk,82.53.Kp 



I. INTRODUCTION 

The simplest description of Bose-Einstein condensa- 
tion (BEC) considers an ensemble of indistinguishable 
particles all in the same quantum state P, Such 
a picture leads to a mean-field representation of the 
many-body wave function as a direct product of the 
single-particle wave functions. This framework has 
been successful in describing many aspects of BEC 
physics. Mean-field approaches neglect two-body cor- 
relations which arise due to the long range interaction 
between two particles. In a microscopic theory, pair cor- 
relations can be incorporated explicitly by expanding the 
correlation functions into cumulants [3, 3 . Alternatively, 
one can also work directly with the correlation functions 
[^,@ . Can one envision directly measurable consequences 
of the pair correlations? We have recently answered this 
question in the affirmative combining our previous 
work on short-pulse photoassociation |8Ml3| with a treat- 
ment of many-body pair correlations [3|-|6[. Here we dis- 
cuss in detail how the choice of experimentally tunable 
parameters determines feasibility of the proposed exper- 
iment. 

Photoassociation assembles two atoms to a molecule 
with laser light [l^. Considerable experimental effort 
has already been devoted to photoassociation with short 
laser pulses. Early experiments were carried out at room 

In the ultracold 
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temperature in mercury vapor 
regime, first experiments aimed at femtosecond photoas- 
sociation of rubidium dimers have led to dissociation of 
molecules created by the trap lasers rather than associa- 
tion [22, [2^ . Recently, evidence of coherent formation of 
molecules in the electronically excited state was provided 
in a pump-probe experiment with chirped femtosecond 
laser pulses [23 - |28| . The proven experimental ability of 
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pulsed photoassociation paves the way for developing a 
pump-probe spectroscopy to study ultracold two-atom 
correlation dynamics Q. 

Photoassociation takes place at the same inter-atomic 
distances where pair correlations are significant. At ul- 
tralow temperatures, photoassociation spectroscopy with 
continuous- wave (cw) lasers is sensitive to the position of 
the nodes of the scattering wave function describing two 
colliding atoms ^ISj . The two-atom scattering wave func- 
tion is closely related to the reduced pair wave function 
characterizing the two-body correlations 0] . The ob- 
served modulation in the spectroscopic features therefore 
serves as a probe of the static inter-atomic correlations. 
If one is to measure two-body correlations dynamically, 
a non-stationary initial state has to be generated. This 
can be achieved by interaction with an external field 
which is fast relative to the ensuing dynamics. There- 
fore the short-pulse photoassociation scenario [1, Q can 
be employed to generate a non-stationary initial state by 
a purrip pulse and then follow its dynamics by a probe 
pulse Q ■ The idea, sketched in Fig. [U is based on impul- 
sive excitation where a selected part of the ground state 
probability density is suddenly removed to the excited 
potential energy surface [l^. As a result the ground 
state phase space density is no longer stationary. The 
void in the ground state probability density is termed 
the dynamical 'hole' 0, US]- The impulsive excitation by 
the pump pulse also transfers momentum to the ground 
state wave function. The 'hole' dynamics can be moni- 
tored by a suitably time-delayed probe pulse, detecting 
the enhancement of ground state density at a specific lo- 
cation by resonantly enhanced multi-photon ionization 
(REMPI) , cf . Fig. [U or by a combination of probe and 
ionization pulse |7| . 

As in any stroboscopic measurement, the pump and 
probe pulses are required to act on a timescale shorter 
than the observed dynamics [l^ [l^. The correlations 
in ultracold atomic gases are due to threshold effects, 
where the kinetic energy of the colliding atoms is smaller 
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or equal to the potential energy. The timescale of the 
correlations is related to the period of the last bound 
level, Tyib w or to the scattering length as, 

T « fias/hk 0. For the example of ^^Rb with a scat- 
tering length of « 100 oo, the correlation timescale is 
typically of the order of 100 ps to 1000 ps. This timescale 
defines an upper limit to the pulse duration. A more 
quantitative estimate of this limit will be derived below 
by comparing the change in kinetic energy to the spectral 
bandwidth of the pulse. On the other hand, the pulse du- 
ration cannot be too short since then the spectral band- 
width would be very broad. This would result in the 
dynamical 'hole' being smeared over a large range of dis- 
tances, hampering precise control of the 'hole' location. 
Therefore the optimal pump pulse is a compromise be- 
tween impulsive excitation and sharp spectroscopic fea- 
tures. 

The insight gained in the photoassociation studies @- 
[Tsj can be used to determine the pulse parameters when 
constructing a spectroscopic probe of the two-atom cor- 
relation dynamics ■ The choice of the pulse parameters 
controls the location and the momentum of the dynam- 
ical 'hole' 0, [13; [13] ■ It can be analyzed by project- 
ing onto the stationary energy cigcnfunctions. From this 
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FIG. 1: (Color online) Pump-probe spectroscopy of the dy- 
namical hole: A photoassociation pump pulse excites popula- 
tion from the ground state to the first excited state, leaving a 
hole in the ground state wave function and transferring mo- 
mentum to it. The ground state wave function (bottom panel) 
is depicted before the pulse (black solid line), at tmax + '24ps, 
i.e. just after the pulse (blue dashed line) and at tmax + 126 ps 
(red dotted line). For red detuned pump pulses, the hole is 
accelerated toward shorter distances which leads to an en- 
hancement of population at short distances. A suitably time- 
delayed REMPI probe pulse can measure this population en- 
hancement. The effect of the probe pulse is modeled by a 
window operator (shown in thin orange lines) where different 
fillings reflect different probe central frequencies. Another op- 
tion is to choose the same central frequencies for pump and 
probe, requiring an additional REMPI pulse, cf. Fig. 1 of 
Ref. m. 



viewpoint, the pump pulse redistributes the initial scat- 
tering state to bound molecular levels and other scatter- 
ing wave functions. The 'hole' location is determined by 
the detuning from the atomic resonance [3, [l3l ■ The mo- 
mentum of the 'hole' can be controlled by proper chirp 
parameters [l3, [3 and by the choice of an attractive or 
repulsive potential of the electronically excited state [l^] , 
i.e. by the pulse detuning. 

The paper is organized as follows. Section |TT] intro- 
duces the theoretical description of the excitation by the 
pump pulse and 'hole' dynamics in a two-state picture 
as well as the model of the probe pulse absorption by a 
window operator. The effect of the pump pulse is studied 
in Section Hill and the observation of the 'hole' dynamics 
is detailed in Section ITVl Section fVl concludes. 



II. MODEL 

We consider two scattering atoms interacting via the 
electronic ground potential. The center of mass motion is 
integrated analytically and omitted from this description. 
For sufficiently low temperatures only s-wave encounters 
contribute. The ground surface supports both bound lev- 
els and scattering states. A stationary initial state is con- 
sidered. For a BEG this state can be well approximated 
by a scattering state [3, [3] with momentum determined 
by the chemical potential. For a thermal state the initial 
state is a thermal mixture of scattering states In 
both cases the two body correlation decay at long inter- 
atomic distance due to the vanishing interaction. This 
means that for both cases the simulation can be based 
on an initial scattering state with a finite cutoff in in- 
ternuclear distance. A photoassociation laser couples the 
electronic ground state to an electronically excited state 
where the potential can be either repulsive or attractive, 
cf. Fig.m 



A. Excitation by a pump pulse and dynamics 

The time-dependent Schrodinger equation is solved for 
two electronic states, 



^h^^\m)=^\m) and (^I^W) - (j^l^.'^^j 



with the Chebychev propagator [30[. The Hamiltonian 
in the rotating-wave approximation reads 



H 



fiEoS{t) 



aEoS*it) t + T/exc(R)-?lAL, 



(1) 



where T denotes the kinetic energy, Vg/^^^ the ground 
and excited state potential energy curves, and ft the tran- 
sition dipole operator. Since the excitation occurs at long 
range, the i?-depcndcncc of p, can be neglected. 
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The pulse parameters appearing in Eq. ([TJ are Eq cor- 
responding to the maximum field amplitude, a slowly 
varying shape of the laser pulse S{t), taken to be Gaus- 
sian, and the detuning Al- The latter is given in terms 
of the atomic resonance line and the central laser fre- 
quency, Al = ujat ^ The pulse energy is calculated 
from the maximum field amplitude and transform-limited 
full-width half-maximum (FWHM) r. 



with 



eocVTTTr- 
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with Eq the electric constant, c the speed of light and 
r the radius of the laser beam [r = 100/xm is assumed 
throughout this work). 

For blue detuning (Al > 0) excitation into the upper- 
most repulsive potential correlating to the 5s -I- 5p3/2 
asymptote, the 0+ (5s -I- 5^3/2) state, is assumed. For red 
detuning (Al < 0), V^xc is taken to be the 0" (5s + 5p3/2) 
state, well-known for its purely long-range attractive 
well. In both cases, corresponds to the lowest triplet 
state, a'^E+(5s -I- 5s). The potentials at short range are 



found in [31|, |32| . They are connected to analytical long- 
range dispersion potentials (^/R + Cq/R -I-Cs/R with 
coefficients taken from [H, [SJ] . We restrict our model 
to this two-state description, considering only the low- 
est triplet state as electronic ground state. In principle, 
a laser pulse excites scattering amplitude from both the 
lowest triplet state and the singlet ground state into all 
states correlating with the 5s-|-5j53/2 asymptote. The res- 
onance condition might differ for singlet and triplet, giv- 
ing rise to different probe signals for singlet and triplet. 
As already pointed out in Ref. , this effect may be used 
to map out the pair correlation functions for the singlet 
and triplet components. 

The Hamiltonian is represented on a Fourier grid with 
variable step size |35l - [37| . Our grid extends to about 
-Rmax = 28000 ao. This ensures that no reflection of the 
wave packet at the outer boundary of the grid occurs: 
Within 10 ns, the fastest wave packet components occur- 
ring in our simulations reach about R = 800 ao. A single 
scattering state with a scattering energy corresponding 
to 20 /xK is considered as initial state in the calculations 
below. The effects of averaging over many thermally pop- 
ulated scattering states and the influence of temperature 
have been discussed in our previous work Q. Here we 
focus on how the pulse parameters need to be chosen in 
order to realize a pump-probe spectroscopy that yields a 
sufficiently high signal to be experimentally feasible. 



W(R)=7^(rpi?,,o)2e-2^'W-^A'(R) 



(2) 



The window operator contains the probe pulse parame- 
ters FWHM Tp and maximum field amplitude Ep^Q. The 
probe pulse central frequency determines the difference 
potential. 



l^{R) = Vp{R)-V^{R)-niCp 



(3) 



with V^(R) denoting the potential which is accessed by 
the probe. The i?-dependence of the transition dipole 
moment fXp is neglected. The physical concept of the 
window operator is to collapse the observation process 
which is completed in a time proportional to the probe 
pulse duration Tp to a single instant in time tp. The finite 
width in time which corresponds to a finite width in fre- 
quency Acj (the bandwidth of the pulse) transforms into 
a finite width in coordinate via the resonance condition 
given by the electronic potentials. This collapse of the 
measurement process assumes that the nuclear motion is 
frozen during the observation, i.e. the excitation is im- 
pulsive, and the window operator W is independent of 
the state of the system [s^. 



III. INDUCING THE CORRELATION 
DYNAMICS BY EXCITATION WITH THE 
PUMP PULSE 

The dynamical pair correlations are initiated by an im- 
pulsive photoassociation pulse. The resonance condition, 
A(i?L) = Vcxc(-Rl) - Vg(i?L) - hwL = 0, determines the 
central position of the dynamical 'hole', i.e. the Condon 
radius Rl- The impulsive conditions on the pulse du- 
ration T are determined by this point R^. An estimate 
is based on the requirement that the energy gain due to 
acceleration during the pulse, Aitkin, is smaller than the 
energy bandwidth of the pulse, ftAw = K/t. For a weak 
pulse such that excitation but no Rabi cycling occurs at 
Rl during the pulse, the energy gain is estimated by the 
acceleration due to the difference potential A(i?), at the 
point of resonance, Rl- In a semi-classical picture, the 
gain in kinetic energy, is estimated by integrating over 
the force, 



P = - 



R=Rl 



B. Modelling the absorption of the probe pulse by 
a window operator 

In order to avoid a separate calculation for each pump- 
probe time delay, the total absorption is represented by 
a window operator W(R) [ssl. Isoj. 

AS = -hiuL^N^ ~ -;iwL(i^g(t)|w(R)|Vg(0) , 



which yields the change in momentum, AP 



Ml 

dR\R=RL 



Evaluating Ai?kin < fiAw leads to an es- 



timate for the upper limit of the pulse duration t. 
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FIG. 2: Upper and lower bounds on the transform-limited 
pulse duration r according to Eqs. © and 



For excitation in the asymptotic region where the differ- 
ence potential can be approximated to leading order by 
—Cz/R^, this becomes 



'asy 



< 



9Ci 



(5) 



For the 0;(P3/2 



state of ^^Rb one obtains r < 85 ps, 
300 ps, and 440 ps for Rl = 50 ao, 78 ao, and 90 ao. A 
lower limit to the duration of the pulse r can be estimated 
due to the requirement that the pulse does not excite the 
atomic transition, i.e. only pairs are excited. This means 
Mw<|14xc(i?L)|«|C3/i?i|or 



> 



hRl 



(6) 



For the 0~(P3/2) state of °'Rb one obtains r > 385 fs, 
1.3 ps, and 2ps for Rl = 50 ao, 78 ao, and 90 ao. The 
bounds on t are plotted in Fig. [21 They indicate how to 
determine the optimal transform-limited pulse duration 
of the photoassociation pulse as a compromise between 
impulsive excitation and sharp spectroscopic features. 
Note that the upper limit scales with while the 

lower limit scales with R^ . The different scaling with Rl 
implies that there is a maximum radius i?2 for which a 
'hole' can be drilled in the pair correlation function. How- 



ever, this maximum radius amounts to R*, 



160000 ao 



or i?2 ~ 8.5 ^m with a corresponding r of 10 ms, so it 
does not impose any practical limitations. The region of 
interest for R^ can be read off Fig.[T] The minimum i?^ is 
due to the decrease in scattering amplitude with shorter 
intcrnuclcar distance. The maximum R^ is determined 
in terms of the pulses that can feasibly be produced in 
an experiment ~ larger Rl require longer delay stages 
between pump and probe pulses and larger transform- 
limited pulse durations, cf. lower panel of Fig. [21 or 
smaller bandwidths, respectively. In the following cal- 
culations, we choose a transform-limited pulse duration 
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FIG. 3: (Color online) The effect of the pump pulse: Pop- 
ulation transfer to the excited state (top), formation of 
ground state molecules (middle) and momentum kick (bot- 
tom). Shown are the excited state population (top), the pop- 
ulation in bound levels of the electronic ground state (middle) 
and the ground state momentum (bottom) after the pump 
pulse is over as a function of the pulse energy. Different detun- 
ings are compared with corresponding Condon radii of 90 ao 
(circles), 78 ao (triangles), and 50 ao (diamonds). 



of the photoassociation pulse of t = 10 ps corresponding 
to an energy bandwidth of 1.5 cm^^. 

The photoassociation pump pulse transfers population 
to the electronically excited state, cf. upper panel of 
Fig. [3l For small and medium detuning, Rabi oscilla- 
tions are observed in the population transfer to the ex- 
cited state, Pexc(ifinai)- The difference between the data 
for Al = — 2.7cm~^ and all other curves is explained by 
excitation of the atomic resonance: the amplitude of a 
10 ps pulse is stiU sufficiently large at zero detuning to 
excite population. The pump pulse does not only trans- 
fer population to the excited state, it has also an effect 
on the ground state wave function, creating the hole and 
transferring momentum, cf. lower panel of Fig. [31 The 
hole corresponds to population of both bound levels and 
higher energy scattering states, i.e. to the formation of 
ground state molecules and hot atom pairs. The popula- 
tion of bound levels is shown in the middle panel of Fig. [31 
For small and medium detuning, molecules are formed 
mostly in the last bound level, w = 40 (between 90% and 
99%). The small amount of initial population near the 
Condon radius is the limiting factor for the population of 
bound levels at large detuning. In that case, molecules 
in the last three levels are formed, with 55% to 65% in 
the level w = 39. These levels can be thought of as mak- 
ing up the hole, i.e. they determine the hole dynamics. 
The bound population decreases at large pulse energies 
for small detuning. This is due to power broadening, i.e. 
population residing at the 'right' distance is transferred 
to continuum states rather than to bound levels. 

The momentum expectation value of the ground state 
wave function after the pulse is shown in the bottom 
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FIG. 4: (Color online) Formation of ground state molecules 
due to the pump pulse. Shown is the projection of the final- 
time wave packet onto the bound ground state vibrational 
levels for two red detunings. 



panel of Fig. [31 For blue detuning, the hole is accelerated 
toward larger distances, while for red detuning dynamics 
toward shorter distances occur (as indicated in Fig. [T|). 
Population enhancement at short distance which is mea- 
sured by the probe pulse can therefore only be expected 
for red detuning. 

While some bound levels arc populated, population 
transfer occurs also into continuum states, and the over- 
all energy of the wave function is changed due to the 
interaction with the pump pulse. For red detuning, an 
increase in energy is always observed, while for blue de- 
tuning a decrease in energy occurs for medium and large 
detuning at large pulse energies. The projection of the 
wave packet onto the bound ground state vibrational lev- 
els after the pulse is shown in Fig. S] for red detuning 
(Al = — 4cm~^). The vibrational distribution will be 
used to interpret the time-dependent signals for weak- 
field excitation. Dynamical signatures of the levels with 
binding energy of Icm^^ and less are expected. The 
corresponding vibrational periods are between 59 ps and 
3 ns. 



IV. OBSERVATION OF THE CORRELATION 
DYNAMICS BY THE PROBE PULSE 

The probe pulse measures the two-atom amplitude at 
a certain range of distances. The dynamics arc unrav- 
eled by varying the delay time between pump and probe 
pulse. Two-atom correlations may be detected by prob- 
ing the molecular part of the ground state wave function. 
Detection of weakly bound ground state molecules has 
been developed in the context of photoassociation exper- 
iments using cw excitation [40| . It is based on resonantly 
enhanced multi-photon ionization. In the following, two 
different detection schemes are explored. 

(i) The dynamics of the hole can be probed when it 



arrives at short internuclear distance in the inner region 
of the potential. This implies the two-photon ionization 
scheme that has been utilized before in photoassocia- 
tion experiments to detect molecules in the lowest triplet 
state [3^ ■ The resonant enhancement is provided by the 
■^Sg (5s -I- 4d) state, cf. Fig. [U i.e. the probe pulse has a 
central frequency that is different from the pump pulse. 
Assuming ionization from the ■^I]+(5s -f- 4(i) intermediate 
state to be saturated, the detection is determined solely 
by the Franck-Condon factors between the lowest triplet 
and the intermediate states. This is reflected by the cor- 
responding difference potential that enters the window 
operator, Eq. 

(ii) Alternatively, the hole can be probed at the posi- 
tion of its creation. A possible REMPI scheme consists 
of overlapping a probe pulse that has the same central 
frequency as the pump pulse with another pulse ionizing 
from the 0" (5s -|~ 5^3/2) excited state. Such an ionization 
scheme has been utilized before in femtosecond photoas- 
sociation experiments to detect excited state molecules 
p^ ]. Assuming the detection to be determined by the 
probe pulse, the difference potential of the 0" (5s -I- 5^3/2) 
and the lowest triplet states enters the window operator. 

In the second detection scheme, one needs to ensure 
that the excited state wave packet that is created by the 
pump pulse does not interfere with the density that is ex- 
cited by the probe pulse. This can achieved by sending 
another ionization pulse simultaneously with the pump 
pulse such that any initial excited state density is imme- 
diately eliminated from the sample. The overall pump- 
probe scheme consists therefore of two time-delayed pulse 
pairs, pump -I- ionization pulse followed by probe -I- ion- 
ization pulse. Since all pulses just address wave packet 
density in a certain range of internuclear distances, no 
phase relation of the four pulses is required. 



Probing the hole at short and intermediate 
distances 



Since the excitation by the pump pulse leads to a re- 
distribution of the ground state probability, the hole con- 
tains a part corresponding to bound molecules and a part 
corresponding to hot atoms. Overall, the hole is acceler- 
ated toward short distances for red detuning of the pump 
pulse. The signal is therefore expected to consist of two 
parts: At early times, a large peak reflects the arrival 
of both molecular and hot atomic density in the probe 
window. After the first reflection at the inner turning 
point of the potential, the hot atomic density leaves the 
short-range distances. Only the molecular part continues 
to move within the range of the potential. Small oscil- 
lations refiect the motion of the molecular part at later 
observation times. 

This behavior is indeed observed in the time-dependent 
expectation value of the window operator modeling the 
probe pulse absorption in Figs. [5] and [T] (W(i)) is shown 
as a function of the pump-probe delay for probe wave- 
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FIG. 5: (Color online) Probing the dynamical hole near the 
inner turning point: The time-dependent expectation value of 
the window operator is plotted for different energies and de- 
tunings of the pump pulse and a probe wavelength of 680 nm 
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FIG. 6: (Color online) Spectrum of the time-dependent expec- 
tation value of the window operator (middle panel of Fig. [S] 
i.e. Al = — 4cm~'^ and a probe window around R — 10 ao). 
The black arrows indicate the position of ground state vibra- 
tional levels, the grey arrows some of the beat frequencies 
between different vibrational levels. 



lengths of 680 nm with 20 nm bandwidth and of 518.5 nm 
with 1 nm bandwidth, respectively. The spectra of these 
signals are obtained by filter-diagonalization [4l|, |42j, a 
method allowing to accurately extract frequencies from 
just a few oscillation periods. The spectra of the probe 
signals of the middle panels of Figs. [5] and [7] are shown 
in Figs. [5] and [HI 

For the probe measurement at = 10 ag (wp = 
680 nm), cf. Fig. [5l an enhancement of the signal by 
a factor up to 10 is observed at short delay times. The 
enhancement is larger for smaller pump pulse detuning 
since more population is excited creating a larger hole, 
cf. Fig. [3l It is due to the accelerated population which 
contains both bound and continuum parts. At longer de- 



lay times, only the bound part of the wave packet returns 
to i? = 10 ap. This part of the dynamics is characterized 
by a large oscillation with overtones. It can be analyzed 
in terms of the vibrational periods of the molecular levels 
making up the bound part of the wave packet. The oscil- 
lation period corresponds to the vibrational period of the 
last, respectively last but one, level, cf. Fig. |4l The fre- 
quency of the long-term oscillations increases with pump 
pulse detuning. This may be rationalized in terms of the 
vibrational distributions, cf. Fig. 

The overtones represent a beating between the differ- 
ent vibrational levels making up the hole. This becomes 
evident in Fig. |6] where the spectrum of the probe signal 
is plotted for the detuning Al = — 4cm~^ and differ- 
ent energies of the pump pulse. The solid arrows in- 
dicate the spectral positions of the ground state vibra- 
tional levels, except for the last bound level which is too 
close to zero to be visible {E]^'^^^ = 8.4 x 10~''cm~^ or 
25 MHz). Of the remaining levels, the second to last 
level at i^Jj^nV^ = 2-1 x lO-^cm^i draws the highest 
peak for weak pump pulses in accordance with the vi- 
brational distribution, cf. lower panel of Fig. H) In the 
weak-field regime, upper panel of Fig. [6l the peaks in 
the spectrum appear at the same positions for different 
pump pulse energies. Increasing the pulse energy leads to 
more population in the lower lying vibrational levels and 
hence additional features, i.e. beat frequencies, appear in 
the spectrum. The strong-field and weak-field regime are 
compared in the lower panel of Fig. [51 At a pump pulse 
energy of 23.5 nJ, Rabi cycling occurs and significantly 
more momentum is transferred to the ground state wave 
packet than for 3.8 nJ, cf. lower panel of Fig. |31 This re- 
sults in a shift of the spectral peak positions, cf. the solid 
blue and dashed green lines in the lower panel of Fig. |3l 
For strong field, so much energy is pumped into the sys- 
tem that the dynamics is not easily unraveled in terms 
of a decomposition into field-free vibrational periods. 

The overall probe signal is much larger for a measure- 
ment at R ~ 24 ao {ujp = 518.5 nm ) than for one at 
R = 10 ao (wp ~ 680 nm). However, the relative en- 
hancement is significantly smaller, cf. Fig. [T] Since the 
measurement occurs at larger distances, there is ground 
state amplitude initially within the probe window, cf. 
Fig. [TJ leading to a strong background. The pump pulse 
cycles population which shows up as a dip in the probe 
signal around zero delay, cf. Fig. [T] Due to the shape of 
the potentials, the probe window is much broader for 
LUp = 518.5 nm than for cjp = 680 nm for the same 
probe spectral bandwidth. The oscillations and over- 
tones at large delays are therefore less clearly resolved for 
ujp = 518.5 nm. Similarly to Fig.[5j larger pump pulse de- 
tuning leads to faster oscillations and larger pump pulse 
energy to a stronger signal in Fig. [T] 

The probe spectra for a measurement at i? = 24 ao 
{u!p = 518.5 nm ) and a pump pulse detuning of Al = 
— 4cm~^ are analyzed in Fig. (5] The signal background 
leads to a non-zero offset of the spectral baseline. The 
positions of the peaks are again compared to those of the 
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FIG. 7: (Color online) Probing the dynamical hole at interme- 
diate distances: The time-dependent expectation value of the 
window operator is plotted for different energies and detun- 
ings of the pump pulse and a probe wavelength of 518.5 nm 
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FIG. 8: (Color online) Spectrum of the time-dependent expec- 
tation value of the window operator (middle panel of Fig. [T] 
i.e. Al = — 4cm"~^ and a probe window around R — 24ao). 
The solid arrows indicate the position of ground state vibra- 
tional levels. 



vibrational levels, indicated by black arrows. As expected 
they agree very well, with the beating between different 
vibrational levels leading to a splitting of the peaks. For a 
strong pump pulse, similarly to Fig.[n]a shift of the peak 
positions is observed compared to smaller pump pulse 
energies. The similarities in Figs. IH] and [51 i.e. peak 
positions and shifts, reflect the underlying wave packet 
dynamics which is completely determined by the (identi- 
cal) pump pulse excitation. The differences such as the 
spectral baseline and specific peak shapes are attributed 
to the differing ways of measuring the dynamics. 

The third window operator shown in Fig. [T] that corre- 
sponds to a probe wavelength of 516.5 nm and to a probe 
bandwidth of 1 nm does not show any time-dependence 
(data not shown). The population over the broad range 
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FIG. 9: (Color online) Probing the dynamical hole at the po- 
sition of its creation: The time-dependent expectation value 
of the window operator is shown for different pump pulse 
energies and a probe wavelength of 780.5 nm (window local- 
ized around i? = 90ao for Al ~ — 2.7cm~^, i? = 78ao for 
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of distances covered by this window does not sufficiently 
vary. 

In conclusion, probing the correlation dynamics by ion- 
ization via the '^E+(5s Ad) state allows for identifying 
the positions of the last bound levels of the ground state 
potential, provided time delays of a few nanoseconds can 
be realized. Probing the dynamics at very short distance, 
i? = 10 ao, yields a signal that clearly oscillates as a func- 
tion of the pump-probe time delay. The strength of the 
probe signal is, however, somewhat discouraging. While 
a larger signal can be achieved by probing the dynam- 
ics at intermediate distance, R = 24 ao, the contrast of 
the oscillations is diminished due to the broader probe 
window. 



B. Probing the hole at the position of its creation 

A similar time-dependence of the probe signal is ex- 
pected for probing the correlation dynamics at the posi- 
tion where the hole is created and at short internuclear 
distances, i.e. a large enhancement at short times due 
to both the hot atomic and the molecular components 
of the hole and smaller oscillations at longer times that 
can be attributed to the molecular components. Addi- 
tionally, when probing at the position of the hole cre- 
ation and with pump and probe pulses that overlap in 
time, around zero delay, a strong dip of the probe sig- 
nal should indicate transient population transfer to the 
excited state. These expectations are confirmed by in- 
specting the probe signals shown in Fig. [31 Moreover, a 
large signal strength and a strong relative enhancement 
or, respectively, depiction are observed. 

The motion of the ground state density toward the 
inner turning point followed by the refilling of the hole 
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FIG. 10: (Color online) Spectrum of the time-dependent ex- 
pectation value of the window operator for Al = — 2.7cm"^ 
and a probe window around ii — 90 ao, Al = — 4cm^^ and 
a probe window around ii — 78 ao (upper panel) and for 
Al ~ — 14cm~^ and a probe window around ii = 49ao (lower 
panel). The black arrows indicate the position of ground state 
vibrational levels, the gray arrows some of the beat frequen- 
cies between different vibrational levels. 



is reflected in the peak of the signal which occurs, de- 
pending on the pump pulse energy, between 840 ps and 
1100 ps for Al = -2.7 cm" \ between 480 ps and 540 ps 
for Al = — 4.0cm~^, and between lOOps and llOps for 
Al = — 14cm~^. The refilling of the hole corresponds 
to the 'recovery of the bleach' known from molecular 
pump-probe spectroscopy. Unlike those examples from 
gas-phase or condensed-phase dynamics, the signal in our 
case is not caused by vibronic motion of molecules but 
rather by the two-body correlations between atoms that 
are present in ultracold gases. 

The highest pump pulse energy always yields the 
fastest refilling. This is attributed to the momentum 
transfer which increases with pulse energy, cf. Fig. |31 
After the initial refilling, another depletion is observed 
followed by a slow recovery with small oscillatory mod- 
ulations of the signal. As in Section IIV Al a larger de- 
tuning of the pump pulse leads to a higher frequency 
of these oscillations. The amplitudes of the oscillations 
can be enhanced by increasing the pump pulse energy, in 
particular for large detunings. 

The strength of the probe signal without any pump 
pulse is given by the probability density amplitude at the 
Condon radius. For larger detuning, the Condon radius 
and hence the amplitude decrease, cf. Fig. [TJ This deter- 
mines the strength of the background signal from which 
the 'recovery of the bleach' is measured. It is decreased 
by an order of magnitude when increasing the pump pulse 
detuning from Al = — 4.0cm~^ to Al = — 14cm"^ (cf. 
initial values in Fig.[S]). 

It might seem a little surprising at first sight, that an 



enhancement of the signal above its background value is 
observed for strong field and large detuning, cf. blue solid 
curve in the lowest panel of Fig. |9l This implies that, at 
certain times, corresponding to pump-probe delays of e.g. 
610 ps or 2300 ps, more probability density resides within 
the probe window than is initially there. One would 
rather expect that, after the initial refilling of the hole, 
when the hot atomic component of the ground state wave 
packet has left the short internuclear distances of the 
probe window, the signal could not be completely recov- 
ered anymore. This picture holds, however, only for weak 
pump pulses. For a pump pulse energy of 23 nJ, Rabi 
cycling occurs, cf. upper panel of Fig. [31 It is accompa- 
nied by power broadening. A broadened pulse spectrum 
corresponds to addressing initial ground state density in 
a larger range of internuclear distances. Therefore the 
23nJ-pump pulse at a detuning of Al = — 14cm^^ ex- 
cites not only the peak in the initial ground state wave 
function around 50 ag but also part of the two neighbor- 
ing peaks, cf. Fig. [TJ Since some of this population is 
cycled back to the ground state, it may contribute to 
the bound components of the hole and show up in the 
long-term oscillations of the probe signal. 

The spectra of the probe signals are examined in 
Fig. [TUl The highest peaks occur at the position of the 
last bound levels of the ground state potential that are 
indicated by the black arrows. The peak around zero is 
due to the last level which is not resolved on the shown 
frequency scales (i^b^nd = ^.4 x 10~^cm"^ or 25 MHz). 
It is off the top of the figure in the upper two panels 
for the smaller pump pulse detunings but reduced in the 
lower panel for the large pump pulse detuning. This is 
in accordance with Fig. [4] where for the larger detuning 
a smaller population of the last bound level is observed. 
The smaller peaks occurring at higher frequencies are 
attributed to beat frequencies between different vibra- 
tional levels, some of which are indicated by the gray 
arrows. Increasing the pump pulse energy leads to larger 
peaks and a finer resolution of the spectral features, in 
particular for the larger pump pulse detuning. A larger 
detuning of the pump pulse implies larger components 
in the ground state wave packet with binding energies 
above 0.1 cm^^, cf. Fig. [H This is reflected by peaks oc- 
curring at higher frequencies in the spectrum (note the 
different scales in the upper and lower panels of Fig. I10|) . 
For a very small pump pulse detuning, Al = — 2.7cm"^ 
and probing at fairly long range, R = 90 ao, the dynam- 
ics are rather slow. An observation time of 5 ns is then 
not sufficient to resolve the spectral features of the probe 
signal. 



V. CONCLUSIONS 

A hole in the pair density of an ultracold gas is the 
consequence of sudden unitary population transfer from 
the ground to the excited electronic state. It represents 
a non-stationary state of the two-atom scattering whose 
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evolution in time can be monitored by a suitably chosen 
probe pulse. Such a pump-probe spectroscopy of pair 
correlations can be implemented by slight modification of 
existing experimental setups |24|-|28f. It can be applied 
to Bose-Einstein condensates as well as thermal ultra- 
cold gases. This corresponds to studying the dynamics 
of a pure state vs those of an incoherent ensemble for 
the timescales probed in such an experiment. The main 
effect of thermal averaging is the contribution of higher 
partial waves. This isparticularly prominent in the pres- 
ence of resonances [31 . If a resonance affects different 
scattering channels differently, such as shape resonances 
in either the singlet or the triplet channel of rubidium, 
pump-probe spectroscopy can be used to map out the 
coupled channels pair wavefunction despite the finite life- 
time of the resonance . 

In the present work, we have studied the influence of 
the pulse parameters on the hole dynamics for a pure ini- 
tial state. Creation of a hole can be induced by a short 
pulse which is faster than the timescale of the nuclear dy- 
namics. A practical requirement is that this pulse should 
be transparent at the atomic transition. In that case, 
only pair correlations are dynamically modified. The 
pulse has to be impulsive compared to the timescale of 
relative motion on the electronic ground state. For ul- 
tracold rubidium, the pulse duration should be a few 
picoseconds or, respectively, the spectral bandwidth a 
few wavenumbers. The same process could be applied to 
other species with the timescale modified to maintain the 
impulsive condition. 

The hole represents a non-stationary superposition of 
many scattering states and a few of the last bound levels. 
The shape of the hole in phase space, i.e. its momentum 
and the position, are controlled by the pulse intensity, 
chirp and detuning from the atomic line. As a result 
there is significant experimental flexibility to interrogate 
specific properties of the pair correlation. The timescale 
of the dynamics of the hole can be estimated by compar- 
ing the acceleration due to the difference potential and 
the energy associated to the bandwidth of the pulse. An 
upper bound of the transform-limited pulse duration is 
related to the gain in kinetic energy due to the gradient 
of the difference potential. A lower bound is imposed by 
the restriction not to excite the atomic resonance. For 
all 'hole' positions of interest, this defines a window of 
possible transform-limited pulse durations which for ru- 
bidium arc of the order of a few picoseconds to tens of 
picoseconds. 



To monitor the dynamical evolution of the hole, a sec- 
ond weak pulse, delayed in time from the pump pulse, is 
applied as a probe. Typically the number of pairs is small 
which means a high sensitivity is required. In this study 
a two-photon REMPI scheme is suggested. The probe 
pulse transfers amplitude to an intermediate state which 
then is ionized leading to a high detection efficiency. We 
have modeled the absorption of the probe pulse by time- 
dependent perturbation theory applied to the first step. 
This assumes saturation of the second step in the REMPI 
scheme. The characteristics of the probe pulse duration, 
chirp and central frequency determines the properties of 
the measurement. 

Probing the two-body correlation dynamics yields a 
signal with clear dynamical features. If one is able to 
provide for pump-probe delays of a few nanoseconds, 
spectral features on a scale of less than lcm~^ can be 
resolved. Since pump and probe just address probability 
density at certain internuclear distances, no phase rela- 
tion between the pulses is required, and long time delays 
between the pump and probe pulse can be realized by 
optical delay stages. 

The full power of coherent control can be employed to 
modify the hole. The simplest modification makes use 
of chirped pulses. Chirping the pump pulse changes the 
shape of the hole, while chirping the probe pulse gener- 
ates a measurement of position and momentum in phase 
space with accuracy limited by the uncertainty relation. 

This pump-probe spectroscopy of the pair correlation 
dynamics can be combined with static external field con- 
trol of the initial pair density. Specifically, tuning a mag- 
netic field close to a Feshbach resonance enhances the pair 
density at short and intermediate distances [1^. Obvi- 
ously, this will lead to a larger probe signal strength. It 
remains an interesting open problem to study the effect 
of the Feshbach resonance on the dynamical properties. 
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